###Exercicio 4.2: Cervejas###
cervejas <-c("chope","lata","garrafa","chope","garrafa", "garrafa","lata","lata","nenhuma","lata","garrafa","garrafa", "garrafa","lata","lata","lata","garrafa","lata","chope","nenhuma", "garrafa","garrafa","garrafa","chope","garrafa","garrafa","chope","garrafa","lata","lata")
#Grafico de barras
barplot(table(cervejas))
#Grafico dotplot
dotchart(table(cervejas))

###Exercicio 4.3###

#1.Histograma do dap dos fustes dos caixeitais
caixeta=read.table("caixeta.csv",header=T, sep=",")
str(caixeta)
summary(caixeta)
class(caixeta)
hist(caixeta$cap/pi)

#2.Histogramas da altura das arvores para diferentes caixatais ("local")
par(mfrow=c(3,1))
hist(caixeta$h[caixeta$local=="chauas"])
hist(caixeta$h[caixeta$local=="jureia"])
hist(caixeta$h[caixeta$local=="retiro"])

#3.Ha diferencas entre as estruturas (distribuicao de tamanhos) dos caixatais?
#As medias e distribuicoes diferem nao apenas entre os tamanhos, mas tambem entre a frequencia de ocorrencia em cada local.

### Exercicio 4.4 ###
e.grandis <- read.table("egrandis.csv",header=T,sep=";")
summary(e.grandis)
class(e.grandis)
head(e.grandis)

#1. Utilize o grafico boxplot para analisar o DAP de rvores de E. grandis em funcao das variaveis regiao (regiao) e rotacao (rotacao). 

par(mfrow=c(2,1))
boxplot(dap~regiao,data=e.grandis)
boxplot(dap~rotacao,data=e.grandis)

#2.Avalie a normalidade da altura do conjunto total de rvores com um grafico quantil-quantil contra a distribuicao normal. 
par(mfrow=c(1,1))
qqnorm(e.grandis$ht)
qqline(e.grandis$ht)
##A varivel altura nao segue uma distribuicao Normal.

### Exercicio 4.5 ###
#1.Analise a relacao dap-altura ('dap' e 'h') em funcao do caixetal (local) com a funcao plot, mas somente para as arvores de caixeta (Tabebuia cassinoides). 

caixeta<-read.table("caixeta.csv",header=T, sep=",")
summary(caixeta)
class(caixeta)
head(caixeta)
caixeta1<-caixeta
caixeta1
caixeta1$dap<-caixeta1$cap/pi
caixeta1
head(caixeta1)
nome.spp<-unique(caixeta1$especie)
nome.spp
caixeta1[caixeta1$especie==nome.spp[3],]
par(mfrow=c(3,1))
plot(caixeta1$h~caixeta1$dap,data=caixeta1,subset=local=="chauas")
plot(caixeta1$h~caixeta1$dap,data=caixeta1,subset=local=="jureia")
plot(caixeta1$h~caixeta1$dap,data=caixeta1,subset=local=="retiro")

#2.Verifique a linearidade com a funcao scatter.smooth
par(mfrow=c(1,1))
?scatter.smooth
scatter.smooth(caixeta1$h~caixeta1$dap,span=1/2)

#3. pacote lattice
library(lattice)
xyplot(caixeta1$h~caixeta1$dap|local, data= caixeta)

